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We study the two-point correlation function in the model of branched polymers 
and its relation to the critical behaviour of the model. We show that the correlation 
function has a universal scaling form in the generic phase with the only scale given 
by the size of the polymer. We show that the origin of the singularity of the free 
energy at the critical point is different from that in the standard statistical models. 
The transition is related to the change of the dimensionality of the system. 



The notion of correlation length is of a paramount importance in discrete 
field theoretical models. To define the continuum limit one requires correla- 
tion length to be infinite. This can be achieved by tuning the system to the 
continuous phase transition. 

Consider as an example a spin system. The correlation length can be defined 
by the behaviour of the connected correlation function : 



where the averaging goes over spin configurations on the lattice of size N. The 
delta function selects pairs at a distance r. This function has a well defined 
thermodynamic (large N) limit. In this limit near the transition point G c (r) 
behaves as 
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Abstract 
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G c (r,N) = (-£%-, V (^ - (s)N){*i ~ (s)n)) 
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Gc(r) - ^) (2) 

where £ = £{(3) is the correlation length which depends on a relevant coupling 
constant (3 and u(x) is a quickly decreasing universal function, typically an 
exponential. In the standard theory of continuous phase transitions the be- 
haviour of £ near the critical point is characterised by the scaling exponent v: 

m - \m~ v (3) 

where A(3 = (3 — (3 C is the deviation from the critical value (3 C . Integrating the 
connected correlation function (1) over r one obtains magnetic susceptibility : 



STr-f^ (M 2 )-(M) 2 

E G d r ) = N = X (4) 

where M = J2i s« is the total magnetization. The susceptibility is a second 
derivative of the free energy density f((3, h) with respect to the magnetic field 
h and has a singularity with a critical exponent 7 : 



x(/?) = ag/(/?,fc)~|A/?r (5) 

Comparing with (2) one obtains 

X-r^HA/^- 2 ^ (6) 
which leads to the Fisher scaling relation : 

7 = (2 -V> (7) 

Similar reasoning can be applied to the energy-energy correlation function, 
which gives specific heat cy when integrated over r. In this case one obtains 
the scaling relation a = 2 — vd for the critical exponent of cy = dpf ~ |A/9| _a . 

In other words the singularity of the free energy is directly related to the di- 
vergence of the correlation length. This picture is closely tied with the renor- 
malization group analysis where infinite correlation length implies the scale 
invariance and hence the existence of a critical fixed point. At the critical 
point the free energy scales uniformly in relevant couplings and this yields 
singularity of the free energy. 

In this paper we discuss a different source of the singularity of the free energy. 
We investigate a model with a random geometry. In such models the lattice is 
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dynamical. There is an ensemble of lattices and one sums over all lattices in the 
ensemble. Geometry is defined by providing the distance definition between 
nodes. 



The correlation function is defined in the same way as in the equation (1). 
Now the average is taken over the ensemble of lattices with size N, say with 
N nodes. If there were fields on lattices one should additionally average over 
them, but here we consider the simplest model without the field dressing. One 
should note here that the correlation function in this case is not just a two- 
point function but rather a global correlator since the distance in the delta 
function depends on the whole geometry. Contrary to the quenched geometry 
models the delta function cannot be pulled outside the average brackets in 



In the following we consider the branched polymer model [1-3] . The partition 
function in this model is given as a weighted sum over an ensemble of trees. 
Trees are weighted by one-vertex branching weights. The partition function 
for the ensemble of trees with N vertices is given by : 



where is order of vertex % and s is a one vertex action and C(T) an appro- 
priate symmetry factor of the graph. For some actions s(q) the system has a 
phase transition. As an example we consider here the model with the action : 



The model has been solved in [2]. It has a fourth order phase transition at 
(3 C = 2.47875.... We discuss this particular form of the action to fix atten- 
tion but the presented results concerning the universality and the critical 
behaviour hold also for a much broader class of actions where (9) is satisfied 
only asymptotically for large q and which in particular allow for tuning the 
critical exponents and the order of the transition with the help of one effective 
parameter (3 C [3]. 

The second derivative of the free energy with respect to (3 is finite but has a 
singularity : 



for small A[3. The model has two phases : the tree phase where the Hausdorff 
dimension is two and the bush phase where it is infinite. At the transition 
the Hausdorff dimension is du = (Pc ~ !)/(& — 2). The entropy exponent 7 S 



(1)- 




(8) 



s{q) = log(g) 



(9) 



(10) 
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changes from the generic value 1/2 in the tree phase to 2 — f3 c in the bush 
phase and 7 S = l/d H at the transition [3]. 



The tree phase has highly universal properties. The volume-volume correlation 
function : 



1 N 



n^) = (-E5Ni4 (ii) 



ho 



is for large N given by the scaling formula [4,5]: 



V(r,N) = VNav(x) (12) 
where the universal function is 

v(x) = 2xe~ x2 (13) 
and the scaling argument 



x = (14) 



with d a finite shift which can be neglected in the large N limit. This form 
holds in the whole tree phase. The only dependence on the coupling is in 
the coefficient a = a((3) which rescales the distance r. The average distance 
between points on the branched polymers is 



<r> " = ErV(r,N) = — (15) 

The meaning of the formulae (12) and (15) is that the Hausdorff dimension is 
two as long as a is finite. At the critical point a diverges and the formula (12) 
does not hold anymore. At this point the Hausdorff dimension changes. 

The volume- volume correlation function (11) is not a measure of geometric 
fluctuations. It is called correlation function only by analogy. Its geometric 
meaning is the average number of points at a distance r from a given point. 
In this way the function (11) defines the dimensionality of the system. On 
a quenched geometry the the scaling (12) and (15) would give the canonical 
dimension. 



To speak about geometric fluctuations and correlations one should consider 
a connected correlation function of the type (1) for a geometrical local quan- 
tity. In the branched polymer a good candidate for such a local quantity is 
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any function which depends on a vertex order. In particular one can con- 
sider correlations between q^ or between s((&) [5,6]. The latter is an analog of 
the energy-energy correlation function and integrated over r gives the second 
derivative of the free energy with respect to (5. 

Using similar techniques as for V(r) one can find the large N limit of other 
correlators. We sketch the derivation in the appendix. In particular for the two 
point connected correlator of the vertex operator s(q) = log(g) we obtain : 



1 N 



G c (r,N) = (—J2$\i-j\A s (<li) ~ ( s (q)}N)(s( qj ) - (s(q)) N ) 



(siq)) 2 ^ g{x) r>0, (16) 



N 

where the scaling variable x is again given by (14). The value at r = is 
a iV independent constant which we denote by c = G c (0,N). The universal 
function g has the form : 



g (x) = v "(x) = 4x (2x 2 - 3) exp (-x 2 ) . (17) 
The coefficients a, 5 and c depend on (3. Near the transition they behave as : 



a ~ \A(3\~ 



2(/3 c -2) 



5 ~ I A/31 sfe^i 
(s(g))oo ~ const + I A/3 1 "-2 
c ~ const + \Aj3\0=2 



One should note here that similar correlation function defined in the grand- 
canonical ensemble vanishes identically for r > 0. It is essential that the 
averages, in particular (s(q))N are taken in the ensemble with a fixed number 
of nodes N. The source of the long-range correlations lies in the fact that for 
fixed iV the orders of vertices satisfy the global constraint 

N 

$> = 2(iV-l) (19) 
i=i 



The second derivative of the universal function v"(x) in the formula (17) comes 
about as follows. One can split the correlation function (16) into three terms : 
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G c {r, N) = (—J2 s \i-j\A s (9i) s (Qj) / ^ 

ij 

1 N 



-2(s(q)) N (-J25\ l - 3 \Mri)^ (20) 



1 N 
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Each of them is calculated separately. For r > and large N they correspond 
to three terms in the following sum : 



G c (r, N) » {s(q)) 2 JV(r + 2-) - 2F(r + -) + V(r)} 

= ( S (g))^v / iVa{t;( a ; + 2-^L) - 2v(x + 4=)+ v(x)} (21) 



' (s(q))la5 2 v"(x) 



<N 

where we used the fact that 



(s(q))N = (s(q)) QO + 0(±) (22) 



The interesting point which we would like to emphasize is that the dressing of 
the geodesic line by the operators s(q) at the ends results only in the effective 
shift of the argument r. This finite shift S gives a subleading contribution to 
each of the correlators in the sum (20). However as the leading terms cancel 
out only those subleading terms contribute to the end result. As we show 
below, the singularity of the free energy is directly related to the singularity 
of 5. 

An important issue is the large iV limit. The formula (16) has been calculated 
for large N. In the scaling variable x (14) the correlation function contains a 
part proportional to the delta function and the scaling part. The scaling part 
gives a non-vanishing contribution when the sum over r > is performed for 
N — > oo even though the contributions for individual r vanish : 



]T G c (r, N) » / drG c (r, N) » ( S (g)>^ 2 / dxv"(x). (23) 

r=l 

The second derivative of the free energy density becomes : 
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Fig. 1. Scaling behaviour of the connected correlation function. The points are 
results of the direct computation of G c (r,N) by the Monte-Carlo program. For 
clarity only every third point is plotted. The line is the universal function g(x). The 
argument x = a((3)r /N 1 / 2 . 

oo 

cv = d}f = £ G c (r, N) = c- 25 2 (s(q)) 2 00 . (24) 

r=0 

The first term c comes from N vertices and is a trivial self correlation. The 
second one is an anomalous term coming from correlations of pairs. The in- 
dividual correlations vanish like 1/N. The number of pairs, however, grows 
as N 2 so effectively the correlations contribute in the same order as c. This 
of course is possible only because the correlations are long range ie. they are 
not cut off by any finite scale. In fact it is the latter term which gives the 
singularity of the free energy density. From (18) we see that a goes to infinity 
and 5 to zero at the critical point. The singularity of c is weaker than the 
one in 5 2 . Thus the singularity of the heat capacity cy comes from the co- 
efficient 5 2 . Paradoxically the contribution from S vanishes at the transition. 
The singularity arises from the way the coefficient 5 vanishes and leads to the 
divergence of higher derivatives. 

In the vicinity of the critical point the universal argument x = ar/N 1 / 2 has 
two large factors : N 1 ^ 2 and a that diverge. One has to carefully define the 
limit A(3 — > 0. If one did it naively and sent A/3 to zero at fixed N, the 
average distance (r) (15) would collapse to zero. In fact it does not. At the 
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transition the Hausdorff dimension changes to du = (/3 C — l)/(/3 c — 2) [3] and 
the universal argument of the functions v and g changes to x oc r/N 1 ^ 11 . Also 
the function v and g themselves change the form. Just beyond the transition 
point the geometry collapses. The collapse results from the appearance of an 
additional mass term e~ m ^ r in the correlation function with the non-zero 
mass independent on N. This mass fixes the average distance between points 
on the polymer to 1/m which is independent of N. The fact that the average 
distance (15) does not grow with the lattice size N can be interpreted as an 
infinite Hausdorff dimension. 

The model can be easily generalized to other actions. By a simple modification 
one has a possibility to tune (3 C in the range (2, oo)[2]. For 2 < (3 C < 3 the 
critical behaviour of the coefficients a, 5 and c is still described by the same 
formulae (18). The coefficient a diverges and S gives the the singularity of the 
free energy. The situation changes for (3 C > 3 where both a, 5 approach finite 
constants at the critical point. The singularities of both of them are of the form 
Af]^ c ~ 3 . The singularity of b is inherited by cy ~ Af3^ c ~ 3 . The fact that a is 
finite at the transition means that the argument x of the universal functions 
(12) and (17) preserves the form x = ar/N 1 ^ 2 and the Hausdorff dimension is 
two as in the tree phase. In the limiting case f3 c = 3 the system undergoes a 
strong third order phase transition with the singularity^ /" ~ log |A/5| _1 . 

Let us at the end briefly compare the critical behaviour of the model with the 
standard critical phenomena of statistical models. The general scaling Ansatz 
for a two-point correlation function involves two dimensionless ratios : r/£ 
and r/N 1 ^ 11 . The singularity of the free energy arises when the correlation 
length £ diverges. In our model the scaling function involves only the argument 
r /N l l dH . In other worlds it has no other scale except the size of the system. One 
can say that the model is always critical (£ — > oo). The long-range correlations 
between orders of different vertices are not induced by local interactions but 
rather by the Euler relation (22). The mechanism of the transition is the 
same as in the balls-in-boxes model and relies on appearance of the surplus 
anomaly [7]. There the correlations between number of balls in two boxes 
vanish for large iV (this is in the analogy with the vanishing of G c (r, N) for 
fixed r > and iV — > oo), but this is compensated by the growing number of 
pairs of boxes over which we have to sum over. The anomaly corresponds to 
the singular vertex which changes the geometrical properties of the branched 
polymer and leads to the collapse which is associated with the appearance of 
the non scaling mass in the two point correlation function. 

A similar phase transition between the branched polymer phase and the col- 
lapsed phase is observed in 4d simplicial gravity [10-12]. One can expect the 
same qualitative behaviour of the correlation function (1) for the curvature- 



2 not /" ~ log | A(3\ as written in [2] 
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curvature correlations in 4d simplicial gravity [13]. 
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A 



The correlation functions are constructed by means of the partition function 
of planted, rooted, planar trees. This partition function can be found from the 
following recursive equation [1]: 

Z = e-^F(Z) (A.l) 



where the generating function is 

F(Z) = Y,e~ P ' { - q) Z q ~ 1 ■ ( A -2) 

q=l 

For fj, approaching a critical value jic from above, Z has the following singu- 
larity : 



Z(fi) = Z ((3) - Z 1 ((3)^i + . . . (A.3) 

The critical value fie corresponds in the large N limit to the free-energy 
density of the canonical ensemble [2]. In particular 



(s(q))N = dpf{0) = -dptoW + O(^) (A.4) 

As a function of j3 the free energy jj,c has a singularity at a critical point (3c- 
This singularity appears also in the coefficients Zj of the expansion (A.3) [2,5] : 



Z ~|A/3|%^ (A.5) 
Z 1 ~\AP\^ (A.6) 



The correlation functions in the grand-canonical ensemble can be calculated 
in terms of Z(/i,(3) [4,5] as : 
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EE4)4)W=(i - ^r(tr) 2 ( A - 9 ) 



which then can be transformed by Laplace transform to the canonical ensemble 
with a fixed size N leading to the terms in (20). The results (22) correspond 
to the leading order terms of the Laplace transform calculated by the saddle 
point approximation. 
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